Massively parallel simulations of relativistic fluid dynamics on graphics processing units with CUDA

摘要

Relativistic fluid dynamics is a major component in dynamical simulations ofthe quark-gluon plasma created in relativistic heavy-ion collisions.Simulations of the full three-dimensional dissipative dynamics of thequark-gluon plasma with fluctuating initial conditions are computationallyexpensive and typically require some degree of parallelization. In this paper,we present a GPU implementation of the Kurganov-Tadmor algorithm which solvesthe 3+1d relativistic viscous hydrodynamics equations including the effects ofboth bulk and shear viscosities. We demonstrate that the resulting CUDA-basedGPU code is approximately two orders of magnitude faster than the correspondingserial implementation of the Kurganov-Tadmor algorithm. We validate the codeusing (semi-)analytic tests such as the relativistic shock-tube and Gubserflow.